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Abstract 

In order to solve quantum field theory in a non-perturbative way, Lagrangian lattice 
simulations have been very successful. Here we discuss a recently proposed alter- 
native Hamiltonian lattice formulation - the Monte Carlo Hamiltonian. In order to 
show its working in the case of the scalar 3>f +1 model, we have computed thermo- 
dynamic functions like free energy, average energy, entropy and specific heat. We 
find good agreement between the results from the Monte Carlo Hamiltonian and 
standard Lagrangian lattice computations. However, the Monte Carlo Hamiltonian 
results show less fluctuations under variation of temperature. We address properties 
of the MC Hamiltonian, like a finite temperature window, and scaling properties. 
Also we discuss possible future applications - like quantum chaos in many-body sys- 
tems, the non-perturbative computation of wave functions of elementary particles, 
as well as scattering amplitudes in high energy physics. 



1 Introduction 



Lagrangian lattice field theory has been for the last three decades the most 
successful non-perturbative technique to compute properties of elementary 
particles and to solve QCD. However, some problems have resisted to a solu- 
tion even by this powerful technique. For example, it is difficult to estimate 
wave functions and the spectrum of excited states. Wave functions in con- 
junction with the energy spectrum contain more physical information than 
the energy spectrum alone. Although lattice QCD simulations in the La- 
grangian formulation give good estimates of the hadron masses, one is yet 
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far from a comprehensive understanding of hadrons. Let us take as example 
a new type of hardrons made of gluons, the so-called glueballs. Lattice QCD 
calculations [1] predict the mass of the lightest glueball with quantum number 
J PC = 0++, to be 1650 ± lOOMeV. Experimentally, there are at least two 
candidates: / (1500) and /j(1710). The investigation of the glueball produc- 
tion and decays can certainly provide additional important information for 
experimental determination of a glueball. Therefore, it is important to be able 
to compute the glueball wave function. 

On the other hand, the Hamiltonian formulation seems very suitable to com- 
pute the energy spectrum and wave function of the ground state and excited 
states. Unfortunately, the Hamiltonian formulation applied to many-body sys- 
tems or elementary particle physics has been widely considered to suffer from 
the absence of a solid many-body solution technique. There are a few ex- 
ceptions: In 1-dimension, the density matrix renormalisation group technique 
works very well [2]. In Ref.[3,4] we have made a proposal how to overcome this 
problem, i.e. by constructing an effective Hamiltonian via Monte Carlo - the 
so called Monte Carlo Hamiltonian (MCH). 

In recent years, Monte Carlo methods have been widely used to solve problems 
in quantum physics. For example, with quantum Monte Carlo there has been 
improvement in nuclear shell model calculations [5]. A proposal to solve the 
sign problem in Monte Carlo Greens function method, useful for spin models 
has been made by Sorello [6]. Lee et al. [7] have suggested a method to di- 
agonalize Hamiltonians, via a random search of basis vectors having a large 
overlap with low-energy eigenstates. 

In this article we will not review the mathematical formulation defining the 
Monte Carlo Hamiltonian, but rather discuss the underlying physical idea. 
Also, we will examine some of its characteristic properties, namely the fact 
that the MCH describes physics in a window of energy and temperature. 
We ask - which are the parameters which determine the position and size of 
the window? How can it be widened? Also we examine a scaling property - 
although the MCH is calculated at a certain inverse temperature (3 , it holds 
for all j3 in the respective window. Finally, we discuss avenues of possible 
future application of the MCH - like quantum chaos in many-body systems 
of quantum field theory, wave functions in elementary particle physics, and 
scattering in particle physics. 



2 Monte Carlo Hamiltonian - an effective Hamiltonian 

The use of Hamiltonian methods in many-body physics has a long history. 
In general, the Hamiltonian is an operator - expressed in terms of fields and 
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Fig. 1. Energy spectrum in a low energy window. 



its conjugate momenta. A widely adopted strategy to compute physics from 
such Hamiltonian proceeds via construction of a Fock-space, on which such 
Hamiltonian operates, and to construct a finite dimensional model subspace 
(e.g. Tamm-Dancoff approximation). By computing matrix elements from such 
subspace, one computes a finite dimensional matrix and obtains - via diag- 
onalisation - a spectrum of energies and wave functions. The problem with 
such procedure is that a cut-off in Fock space particle number and in high 
momentum is uncontrolled. A much better idea is the renormalisation group 
idea by Kadanoff and Wilson. They suggest to construct (via blockspin trans- 
formation or intregrating out high momentum components) a Hamiltonian for 
the physics of critical phenomena, valid in the neighborhood of a fixed point 
- which needs much less degrees of freedom than the original Hamiltonian. 

The Monte Carlo Hamiltonian is similar to this in the sense that it aims to 
decribe physics only in a finite range of parameters (i.e. a low energy window). 
In lattice field theory, one constructs a Hamiltonian via the transfer matrix - 
i.e. going to the limit of transition time T — > 0. This leads - for lattice gauge 
theory - to the well known Kogut-Susskind Hamiltonian, which is expressed in 
terms of gauge fields (link variables) and its conjugate momenta. In contrast, 
the MC Hamiltonian is obtained from transition matrix elements at finite 
transition time. As a result, the MC Hamiltonian is obtained and defined in 
terms of its matrix elements - and not in terms of operators. As pointed out 
above, the crucial question is: In which subspace or model space are those 
matrix elements to be taken? Contrary to the Tamm-Dancoff Hamiltonian ex- 
pressed in terms of particle creation and annihilation operators in momentum 
space, the MC Hamiltonian is formulated in a space, which is the analogue 
of position space in quantum mechanics - the so called Bargman space. It is 
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interesting to note that the problem to conserve manifestly local gauge invari- 
ance in lattice gauge - apparently necessary for confinement of quarks - has 
been solved by Wilson, using a position space formulation (e.g. Wilson loops) 
and not a momentum space formulation. It looks like as if the description of 
gauge theory in nature is simpler in position space. 



However, the use of Bargman space is by itself not good enough to solve a 
many-body problem. An idea - essential for the MC Hamiltonian - is to use a 
model subspace created in a random way (Monte Carlo), however, guided by 
a physical principle, which gives a large ponderation to such states which are 
expected to yield "large" (non-negligeable) transition matrix elements. This 
idea has been taken over from Lagrangian lattice field theory, where the solu- 
tion of path integrals - via Monte Carlo with importance sampling - allows to 
compute (estimate) physical observables by summing over field configurations 
$ which give a "large" weight exp[— S(<&)/h] where S is the action. In other 
words, we suggest that Monte Carlo with importance sampling, which allows 
to solve lattice field theory in the Lagrangian formulation, is also the key to 
solve lattice field theory in the Hamiltonian formulation. This generates what 
we call an (importance sampled) stochastic basis in Bargman space, which is 
absolutely essential for the solution of the Hamiltonian many-body problem. 




Fig. 2. Free energy F((3). Comparison of results from Monte Carlo Hamiltonian 
(filled circles) with standard Lagrangian lattice calculations (open circles). 
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3 Model: chain of coupled anharmonic oscillators 



We consider a one-dimensional chain of N s coupled harmonic oscillators, with 
anharmonic perturbation. Its Euclidean action is given by 



In the continuum formulation it corresponds to the scalar model, 



Introducing a space-time lattice with lattice spacing a s and a t , this action 
becomes 



N 3 N t -1 
n=l k=0 



1 / $(x n ,t k+1 ) - ${x n ,t k y 



a t 



+ 



i / $(x n+ i,t k ) - ®(x n ,t k y 



m 



f ^<5> 2 {x n ,t k ) + ^.<5>\x n ,t k ) 



(3) 



The actions given by Eq.(3) and Eq.(l) can be identified by posing <fi = y / al < l ) , 
Q = l/a s , f2 = m, and A/2 = g/4\. 



4 Comparison: Monte Carlo Hamiltonian vs. Lagrangian lattice 
simulation 



Energy spectrum. In Ref. [8] we have tested the method for a chain of harmonic 
oscillators (Klein Gordon model) where the energy spectrum is known exactly. 
We found quite good agreement between the results from MC Hamilton and 
the exact results. The energy spectrum for the chain of anharmonic oscillators 
(scalar model) obtained from the MC Hamiltonian is shown in Fig. [1]. The 
model parameters are: Q = 1, Q = 2, A = 1 (H = 1, k B = 1). Approximation 
parameters are: Lattice size N s = 9, lattice spacing a s = 1, dimension of 
stochastic basis N stoch = 100. The statistical errors (indicated by error bars) 
are reasonably small and are more or less of the same size for n — 1 to n — 40. 
They increase noticeably for n > 40. Also the shape of the curve E n vs. n 
changes beyond n = 40. This signals the occurence of a low-energy window 
stretching from E 1 = 11.27 to E m = 16.52. 
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Fig. 3. Same as Fig. 2, for average energy U{(5). 
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Fig. 4. Same as Fig. 2, for entropy £*(/?)• 

Thermodynamical observables. A solid test of the Monte Carlo Hamiltonian 
method is a comparison with results from standard Lagrangian lattice calcula- 
tions. The information of the energy spectrum enters into such thermodynam- 
ical functions. Thus thermodynamical functions allow an indirect test of the 
energy spectrum. We have chosen to compute the following thermodynamical 
observables: the partition function Z, free energy F, average energy U, specific 
heat C, entropy S and pressure P. Those are defined by 
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Z{/3) = Tr[exp(-/3H)] , 
F(P) = -~]ogZ, 

U(P) = ±Tr [Hexp {-(3H)] = 
S(P) = -(U-F)=k B (3(U-F) , 

T 

dF 



dV 



(4) 



Here k B denotes the Boltzmann constant. The temperature r is related to (3 
via /3 = T/H= l/(k B r). 
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Fig. 5. Same as Fig. 2, for specific heat C{j3). 

Fig. [2] shows the free energy F as a function of (3. The following behavior is 
observed. First, there is good overall agreement in the range 0.5 < (3 < 10. 
However, there is an essential difference between the Lagrangian approach and 
the MC Hamiltonian approach: In the Lagrangian approach, one needs to do 
a new computation for each value of temperature or inverse temperature f3. 
Moreover, the computation is easier in the regime of small (3, i.e. high temper- 
ature. Increasing f3 requires larger lattices, which means that the numerical 
effort and errors increase. In the MC Hamiltonian approach, one compute tran- 
sition matrix elements for a particular fixed value (3q. One extracts a spectrum 
of energy eigenvalues and corresponding wave functions. Suppose we are close 
to the quantum continuum limit. Then the obtained energy spectrum should 
become independent of /3 . Then plugging those energies in the expressions 
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of thermodynamical functions, Eq.(4), allows to obtain those function for any 
value of (3 (asymptotic scaling). Because in practice one can never reach the 
continuum limit nor the infinite volume limit, as a consequence one will not 
be able from a single MC Hamiltonian siumulation to obtain thermodynam- 
ical functions for any (3. Rather one can compute those functions only in a 
bounded domain in (3. In other words there is a finite temperature window. 
This is one essential difference between both approaches. There is another dif- 
ference. While the Lagrangian approach is efficient and precise in the domain 
of small (3 or high temperature, the MC Hamiltonian approach is efficient 
and precise in the domain of high (3 or low temperature. This can be seen in 
Fig. [3], which displays the average energy. With increase of j3, one observes 
that the Lagrangian lattice results fluctuate more than those from the MC 
Hamiltonian. This phenomenon is more evident in Fig. [4] for the entropy and 
even more so in Fig. [5] for the specific heat. From the data of F, U, S and 
C, we estimate the temperature window of the Monte Carlo Hamiltonian to 
range from f3 — 1 to f3 — 10. 

One may ask: On what does the location and size of the finite temperature 
window depend? It certainly depends on the parameters which control the 
approach to the continuum limit and the infinite volume limit, i.e. lattice 
spacing a s and a t , on the lattice size N s and N t . In lattice QCD it will depend 
also on the coupling g s and gt- It will depend on N st0 ch, the size of the stochastic 
basis. Finally, it will depend also on the precision in computing the transition 
amplitudes (statistical error, etc.). From the study of the scalar model, we 
found that the size of the lattice in temperature direction as well the size of 
the stochastic basis are most important. 



5 Outlook: Possible applications of MC Hamiltonian method 

In order to make the MC Hamiltonian a valuable tool in condensed matter 
or elementary particle physics, it is crucial to show that it allows to treat 
gauge theories and eventually fermions. Work on incorporating gauge theo- 
ries is presently under way. If the MC Hamiltonian methods works for such 
theories, there is much interesting physics to be done, where conventional 
non-perturbative methods, in particular Lagrangian lattice field theory, have 
brought little progress. Let us outline a few examples, (i) Energies and wave 
functions of excited states. Here, in our opinion, the Hamiltonian formulation 
is most suitable. Many-body wave functions play a role for hadron structure 
functions in particle physics, electromagnetic form factors in nuclear physics, 
Bose-Einstein condensation in atomic physics, (ii) S-matrix and scattering am- 
plitudes in high energy physics. The non-perturbative calculation of scattering 
amplitudes in high energy physics is an unsolved problem yet. In Ref. [9] Kroger 
has suggested how to compute the S-matrix in a time-dependent manner from 
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spectral information, i.e. energies and wave functions. The MC Hamiltonian 
may provide such spectral information, (iii) Quark-gluon phase transition in 
QCD at finite temperature and finite density. Such phase transition requires 
to treat QCD in the presence of a non-zero chemical potential. In Lagrangian 
Lattice QCD, the treatment of the gauge group SU(3) leads to a complex 
action, which escapes a treatment by Monte Carlo methods [10]. This has 
been a long standing problem. On the other hand, in the Hamiltonian for- 
mulation on the lattice, the presence of the chemical potential does not lead 
to any complex term, (iv) Quantum chaos from level density distributions in 
quantum field theories. It is a widely used strategy to characterize quantum 
chaos, e.g. in nuclear or atomic physics, via level density distributions (Wigner 
for fully chaotic) and Poissonian (for integrabel non-chaotic) systems [11]. To 
obtain good statistics requires the knowledge of a considerable number of en- 
ergy levels. The idea is that such information may be provided by the MC 
Hamiltonian. 



6 Summary 

We have discussed physical ideas underlying the concept of the Monte Carlo 
Hamiltonian. For the case of the scalar model we have computed thermody- 
namical functions and by comparison with standard Lagrangian lattice cal- 
culations we showed that the MC Hamiltonian works. We outlined possible 
future applications of the MC Hamiltonian in elementary particle physics, 
nuclear physics and atomic physics. 
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